The amplitude of sidebranches in the dendritic crystal growth of ammonium chloride 
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We report new measurements for the amplitude of sidebranches in the dendritic crystal growth 
of ammonium chloride from supersaturated aqueous solutions. No simple power law is found to 
describe either the growth of the average sidebranch amplitude or the average sidebranch envelope. 
Instead, the effective power law exponent appears to increase as a function of distance from the 
dendritic tip. We also find that the amplitude of the sidebranches in this system is significantly 
larger than predicted by models of noise-driven sidebranching. 
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I. INTRODUCTION 

Dendritic crystal growth is commonly observed in the 
growth of non-faceted materials from pure melt or super- 
saturated solution. The crystals are characterized by a 
smooth, nearly parabolic tip, with sidebranches emerg- 
ing a short distance behind the tip. Much of the beauty 
of complex crystal structures results from the intricate 
development and subsequent competition of these side- 
branches. For general reviews, see Boettinger et al. [l[ 
and Asta et al. 

Considerable theoretical and experimental effort has 
focused on measuring the properties of the sidebranch 
structure, understanding the origin of the sidebranches 
and looking for scaling laws that might govern their ul- 
timate development. Previous extensive studies of side- 
branch structure have been reported for succinonitrilc 
, ammonium bromide [f| Q , Xenon Q , ammonium 
chloride |9], and succinonitrile and succinonitrile/acetone 
alloys (Tj 

In this paper, we focus on the early and intermedi- 
ate stages of sidebranches in the dendritic crystal growth 
of ammonium chloride at low supersaturation, measure 
their amplitude, and compare to various power-law pre- 
dictions. 



II. THEORY 

The basic theory for steady state diffusion-limited den- 
dritic crystal growth is presented in Langer ll]. The 
crystal is characterized by a smooth, nearly parabolic tip 
of radius p growing at constant speed v. The tip radius 
and velocity are related to the "stability constant" a* by 
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where D is the relevant diffusion constant, and do is the 
capillary length, which is related to the solid-liquid inter- 
face energy. Direct testing of Eq. [T]has proven difficult, 
in part due to challenges in measuring do and D, and in 



part due to the complications introduced by convection 
in most terrestrial experiments. There is also uncertainty 
whether a* is constant, or whether it might vary with in- 
creased undercooling or supersaturation over the range 
of experimental interest (l2|. 



A. Sidebranches 

A short distance behind the tip, sidebranches emerge 
with a characteristic wavelength A that is typically about 
3 ~ 6p. These sidebranches start out approximately uni- 
form, but compete in a complex nonlinear coarsening 
process [!, Q . Larger branches continue to grow, while 
shorter ones stop, or even begin to dissolve back, eventu- 
ally giving rise to structures with a wide range of length 
scales. 

There are two main approaches to theoretically model- 
ing the origin of the sidebranches. One possibility is that 
the tip growth itself is actually slightly oscillatory, and 
the e merging sidebranches reflect that underlying oscil- 
lation |13l4l5j . Another possibility is that sidebranches 
are noise-driven, and their approximate periodicity is the 
result of the selective amplification of that noise [l6[ . 

Directional solidification experiments with pivalic 
acid/coumarin alloys have shown that directly applying 
a pulsing laser to dendritic tips in directional solidifica- 
tion can drive an oscillating tip and produce a sidebranch 
structure with the corresponding wavelength [l7j . Analo- 
gous results have also been obtained in viscous fingering 
experiments [l8|]. These results show that oscillating- 
tip solutions are possible, but leave open the question 
of what happens as the amplitude of the applied noise 
is reduced to zero. For xenon, both noise-driven and 
perturbation-induced sidebranches have been observed 
I- 

For non-axisymmetric needle crystal growth, Brener 
and Temkin [19| predicted that the noise-induced side- 
branch amplitude A(z) is given by 



A(z) = pSo exp 
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where z is the distance back from the tip, w(z) is the 
average width of the dendrite, and So is the dimension- 
less noise amplitude. For ammonium bromide dendrites 
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growing from supersaturated aqueous solution, Gonzalez- 
Cinca et al. [20[ give Sq ~ 6 x 10~ 5 as a conservative 
over-estimate of the noise value. 

More generally, if w(z) scales as a power law with z, 
then Eq. [2] is of the form of a stretched exponential 
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A(z) = pS exp — 
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where a is 0.4 if w(z) ~ z ' , as found by Brener [21|. 
Power law behavior has also been reported for a number 
of integral parameters of dendritic growth, inclu ding the 
sidebranch envelope, contour area, and volume (El I22M2E] . 

Beyond the initial development of the sidebranches, 
a number of approaches have been developed to model 
the interaction and coarsening of sidebranches. Phase 
field models attempt to incorporate the full physics of the 
three-dimensional problem 26, 27] . The addition of noise 
to phase field simulations has also been shown to produce 
sidebranching structures similar to those observed in ex- 
periments (2a - t30l | , though it continues to prove challeng- 
ing to perform the calculations in the small undercooling 
and small anisotropy range appropriate for ammonium 
chloride solution growth [3l|, . 

Other approaches, such as the coupled map lattice 33§ , 
attempt to model the essential features of sidebranch de- 
velopment and competition more directly, and describe 
power-law coarsening behavior seen in sidebranch exper- 
iments dUl]. 



III. EXPERIMENTS 

The experiments were performed with aqueous solu- 
tions of ammonium chloride with approximately 36% 
NH4CI by weight, corresponding to a saturation temper- 
ature of approximately 60° C. The solution was placed 
in a 40 x 10 x 2 mm glass spectrophotometer cell sealed 
with a Teflon stopper held in place by epoxy. The cell 
was mounted in a massive copper block, surrounded by 
an insulated outer aluminum block, and placed on a mi- 
croscope enclosed in an insulating box. The cell was 
mounted horizontally to minimize the effects of convec- 
tion on the growing crystal. For this experiment, we es- 
timate D = 2500^iin 2 /s (35| . similar to the values of 2600 
/xm 2 /s used by Tanaka and Sano [36j], and 2280 pm 2 /s 
used by Sawada et al. (37| . 

The temperature of the outer aluminum block was con- 
trolled to approximately ±0.1°C. The temperature of the 
inner copper block was controlled directly by computer, 
allowing complete programmatic control over the temper- 
ature during the course of a run. The temperature of the 
sample was stable to within approximately ±5 x 10~ 4 °C. 

Images were obtained from a charged coupled de- 
vice (CCD) camera attached to the microscope and ac- 
quired directly into the computer via a Data Translation 
DT3155 frame grabber with a resolution of 640 x 480 
pixels. The ultimate resolution of the images was 0.63 ± 
0.01/im/pixel. 



To obtain a single crystal, the solution was heated to 
dissolve all the NH4CI, stirred to eliminate concentra- 
tion gradients, and then cooled to initiate growth. Many 
crystals would nucleate, but an automated process was 
set up to acquire images and slowly adjust the tempera- 
ture until only a single isolated crystal remained. Once a 
suitably aligned single crystal was obtained, the temper- 
ature was gradually lowered at a rate of 0.75 x 10~ 3 C /s 
to initiate dendritic growth. Once steady state dendritic 
growth was observed, approximately 5300 images were 
acquired at 1-second intervals. 



IV. AVERAGE DENDRITIC CRYSTAL SHAPE 

For each crystal image, the interface position was mea- 
sured as in Ref . [35j . In order to determine the tip posi- 
tion and size, the tip region was fit as in Ref. [38[ by a 
parabola with a small fourth-order correction 
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where z is the perpendicular distance back from the tip, 
x is the dendrite width, and A4 w —0.0036 for this ma- 
terial 38]. Up to distances about 6p back from the tip, 
this fit is fairly robust and does not appear to be unduly 
influenced by the early stages of sidebranching. Simi- 
lar results were found in Ref. [10j for succinonitrile and 
succinonitrile/acetone alloys, though at higher concen- 
trations, the sidebranches tended to encroach upon the 
tip, rendering the fit less robust. 

A typical dendrite tip is shown in Fig. [I] along with the 
fit to Eq.|H and the average width w(z) for all crystals in 
this run. Near the tip, Eq. [4] fits well, but it fails beyond 
about lOp, before the sidebranches become clearly visible. 

We attempted to fit this average shape w(z) to a power 
law. Over the range of distances studied here, namely 
3 < z/p < 80, the best fit power law for the average 
width is given by 
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where c = 0.94±0.01 and (3 = 0.658±0.002. The average 
shape and the power-law fit are shown in Fig. [2j 

As was previously shown in Ref. 38], however, no 
single power law adequately describes the full range of 
sidebranch growth. Near the tip, the shape is nearly 
parabolic, while far from the tip, the largest branches 
grow as independent dendrites. Consequently, the ef- 
fective power law exponent in Eq. [5] increases steadily 
as one moves further from the tip. This is illustrated 
more clearly in the log-log plot of the average width in 
Fig. |31 We did not observe any well-defined region where 
ft = 3/5, as predicted in Ref. 21|. Empirically, we were 



able to get a reasonable fit over this range by considering 
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FIG. 3. Log-log plot of the average width, including fits for a 
FIG. 1. Dendritic crystal of NH4CI. The width of the image single power-law and a superposition of two power laws, 
is approximately 400 fim. The inner curve (blue) is the best- 
fit parabola for the tip. The outer curve (red) is the average 
width of all crystals in this run. The fourth-order fit (Eq.[4]) is 
on the far right in green; the negative fourth-order term causes 
the fit to deviate sharply from the tip before the sidebranches 
become significant. 
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FIG. 2. Average width of a dendritic crystal along with fits 
to a single power law, Eq. [5] and a pair of power laws, Eq. [6] 



a sum of two power laws, 
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where c x = 1.16 ± 0.04, ft = 0.14 ± 0.11, c 2 = 0.46 ± 
0.10, and ft = 0.79 ± 0.03, though the two exponents 
varied with the exact z-range considered. This fit is also 
included in Figs. [5] and 
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A. Sidebranch Amplitude 

We have performed two separate characterizations of 
the sidebranch amplitude. First, we have considered just 
the envelope of active sidebranches, as in Refs. and 
[To| | . A sidebranch was considered "active" if it was larger 
than all other sidebranches on the same side closer to 
the tip. Figure H] shows the active sidebranches for a sin- 
gle crystal and the average envelope for all active side- 
branches. 
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FIG. 4. The active sidebranches for one crystal and the aver- 
age envelope for all images. 

For small z, less than about 15p, detection of individual 
sidebranches is problematic, and the result is quite noisy. 
Over the intermediate range 15 < z/p < 80 covered by 
our data, the best power-law fit for the envelope e(z) is 
given by Eq.[5]with c = 0.63 ±0.01 and /3 = 0.78 ±0.01. 
This fit is shonw on a log-log scale in Fig. [5] Again, 
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FIG. 5. Log- log plot of the measured average envelope for 
all sidebranches, along with the best power-law and stretched 
exponential fits. 

however, the power law is only approximately valid over 
a limited range. The curvature in the measured aver- 
age envelope is clearly evident. Empirically, we were 
able to obtain a better fit by adding a constant term 



Co = 2.3 ± 1.4 to a power law with c = 0.24 ± 0.15 and 
P = 0.97 ± 0.11, though the large uncertainties indicate 
the fit is not well-constrained. Alternatively, we found 
that the average envelope can also be fit by a stretched 
exponential of the form of Eq. G3 with S = 0.44 ±0.11, 
s = 0.24 ± 0.18, and a = 0.22 ± 0.02. (The fit involving 
a power-law plus a constant is nearly indistinguishable 
from the stretched exponential over the data range con- 
sidered here, so it is not included in Fig. [5]) 

As a second measure of the sidebranch amplitude, we 
considered the fluctuations of the width about the aver- 
age. 

Specifically, we first scaled all distances by the tip ra- 
dius p, subtracted the mean shape, and then computed 
the root mean square deviations of the crystal width as 
a function of distance z back from the tip. In order to 
reduce possible issues with measurement noise contami- 
nating the results, we only included the range z > 15p 
in the fits, though the results were insensitive to the pre- 
cise small-z cut-off. At larger distances, when the side- 
branch structure becomes more complex, the rms fluc- 
tuation measurements become less reliable, so we only 
included data with z < 70p. 

The result is shown in Fig. [6] along with the best fits to 
Eq. [5] and [31 The best power-law fit to this data is given 
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FIG. 6. The rms sidebranch amplitude, along with the best-fit 
power-law and stretched exponential curves. 

by Eq.[5]with c = (1.50±0.06) x 10~ 4 and /3 = 2.26±0.01. 
The best stretched exponential fit to Eq. [3] yields So = 
0.0023 ± 0.0015, s = 0.40 ± 0.33, and a = 0.37 ± 0.04. 
The large uncertainties again indicate that the fit is not 
well-constrained by the data. A log-log plot of the data 
and the two fits is shown in Fig. [7] 

Inserting the power-law fit for the average width in 
Fig. [2] into Eq. [2] yields predicted values of s = 0.75 and 
a = 0.49, somewhat larger than the measured values. 
The fit value for the noise amplitude Sq is considerably 
larger than the estimate of 6 x 10~ 5 given in Ref. [2(|, but 
has considerable statistical uncertainty. Over the z-range 
considered here, decreasing So can partially be offset by 
decreasing a slightly and significantly decreasing s, but 
Sq could not be adjusted below about 6 x 10~ 4 without 
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FIG. 7. Log-log plot of the rms sidebranch amplitude, along 
with the best-fit power-law and stretched exponential curves. 



more than doubling the \ 2 value for the fit. 

As a further check on the noise amplitude, we also 
considered the sidcbranches in individual crystal images. 
We modeled the dendrite width for one side of a single 
crystal by 

/ Z \ a ( 27T \ 

w(z) = w(z) + S cxp y-J sin I — z + <f>\ , (7) 

where A is the sidebranch wavelength, is a phase con- 
stant, and all distances have been scaled by p. The 
results for a single crystal are shown in Fig. [8j where 
So = 0.0021 ± 0.0009, s = 0.40 ± 0.15, a = 0.35 ± 0.05, 
and A = 5.12 ± 0.07. These values are consistent with 
those obtained above for all the combined images. 
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FIG. 8. Border of a single crystal along with the best fit for 
a noise-driven sidebranch function. 



V. DISCUSSION 

We have considered three separate characterizations 
of ammonium chloride dendrites, the average width, the 
average envelope, and the rms sidebranch amplitude. In 
contrast to previous studies on xenon Q and succinon- 
itrile alloys (Toj . we do not find simple power- law be- 
havior for any of these quantities. Instead, we find that 
the effective power-law exponent generally increases with 
distance z from the tip. 

We find that a stretched exponential of the form of Eq. 
03 provides a reasonable fit to the sidebranch amplitude 
data, but the value for the noise amplitude So is more 
than an order of magnitude larger than that estimated 
for intrinsic noise [2(| ■ Excluding either the large or small 
z range from the data in Figs. [5] and [7J does not change 
the estimated value for So significantly. The typical in- 
dividual fit in Fig. [8] is also consistent with a large value 
of S . 

The origin of the large noise amplitude is unknown. 
Although the horizontal growth geometry was chosen to 
help minimize convective effects on the sidebranches, it 
is still possible that convection plays an important role. 

It is important to recognize that the range of appli- 
cability of Eqs. [3] and [7] is quite limited. At small z, 
there are two main issues. First, emerging sidebranches 
are potentially masked by measurement noise. Second, 
measurements of the tip itself (and hence all distances 
scaled by p) are also potentially contaminated by early 
sidebranches. These effects are particularly problematic 
near the tip because the concentration gradients and cor- 
responding optical distortions are largest there. Measure- 
ments of the very early sidebranches and tip shape are 
thus inextricably intertwined. At larger z, the increas- 
ing nonlinear competition among sidebranches leads to 
changes in the expected scaling behavior as well as ambi- 
guities in characterizing the structure by a single- valued 
width function w(z). 

These constraints make it more challenging to identify 
whether there are simple underlying scaling laws that 
govern the growth of sidebranches. In future experi- 
ments, we hope to extend the useful range of the data 
and provide tighter constraints on the functional form of 
the sidebranch amplitude. 
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